Effect of strain on hyperfine-induced hole-spin decoherence in quantum dots 
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We theoretically consider the effect of strain on the spin dynamics of a single heavy-hole (HH) 
confined to a self-assembled quantum dot and interacting with the surrounding nuclei via hyperfine 
interaction. Confinement and strain hybridize the HH states, which show an exponential decay for 
a narrowed nuclear spin bath. For different strain configurations within the dot, the dependence of 
the spin decoherence time T2 on external parameters is shifted and the non-monotonic dependence 
of the peak is altered. Application of external strain yields considerable shifts in the dependence 
of T2 on external parameters. We find that external strain affects mostly the effective hyperfine 
coupling strength of the conduction band (CB), indicating that the CB admixture of the hybridized 
HH states plays a crucial role in the sensitivity of T2 on strain. 

PACS numbers: 03.65.Yz, 31.30.Gs, 62.20.-x, 73.21.La 



I. INTRODUCTION 

During the last years, heavy-hole (HH) spins have 
attracted much interest in the field of spintronics and 
spin-based quantum computing. This is because, com- 
pared to the timescales set by the electron spin, very 
long hole spin relaxation times Ti « 1 ms have been 
predictecP^l and confirmed experimentally^. Further- 
more, ensemble hole spin coherence times T2 > 100ns 
have been measurecP. Alongside, the essential ingre- 
dients for processing quantum information successfully, 
hole spin initializatiorP'Sl and coherent control of sin- 
gle hole spin»^"^, have been shown in quantum dots. 
Additionally, methods applicable to overcome decoher- 
ence, e.g. by preparing the nuclear spin bath in a nar- 
rowed state^^lt^, have been introduced. The prolonged 
timescales regarding decoherence are attributed to the 
Ising-like hyperfine coupling of holea^^ due to the p-wave 
symmetry of the Bloch states in the valence band (VB). 
Recently, the hyperfine interaction strength of holes was 
predicted to be approx. 10% of the interaction strength 
of electronsl^. This was confirmed in experi ment s car- 
ried out in self-assembled InAs quantum dot d^^ l ^^l The 
associated hole spin decoherence time T2 was shown 
to depend on external parameters in a non-monotonic 
fashiorPSl. Due to lattice mismatch, the strain profiles 
of InAs/GaAs dots show a compression in the lateral 
plane and a stretching in the vertical directiorpll. The 
associated strain fields are of considerable strength and 
may strongly affect the band hybridization in the dot and 
hence the spin decoherence. For the light hole (LH) and 
HH band, the effect of confinement and anisotropic lat- 
eral strain on band mixing and on the interaction with a 
Gaussian nuclear field distribution via dipole-dipole hy- 
perfine interaction has been considered in Refs. [22 23, 
In the present work, we examine the effect of realistic 
strain distributions on the spin decoherence time T2 of a 
single HH spin confined to a self-assembled InAs quan- 
tum dot interacting with a narrowed nuclear spin bath 
via hyperfine interaction. We follow the procedure out- 
lined in Ref.[2niwith emphasis on the new features coming 



from strain. The emerging band hybridization is strain 
dependent and shows considerable admixtures of the low- 
est conduction band (CB), and the LH and the split-off 
(SO) band of the VB. An effective hyperfine Hamiltonian 
is derived from the hybridized states being, for realistic 
strain configurations, predominantly of Ising form with 
small hole-nuclear-spin fiip-flop terms which cause expo- 
nential spin decoherence. We study the effect of vari- 
ous internal strain configurations and of applied external 
strain on the decoherence rate I/T2 and its dependence 
on external parameters. Applying external strain up to 
the breaking limit of the sample affects the effective hy- 
perfine coupling of the CB admixture much more than 
the coupling of the LH admixture. In contrast to this, 
the changes in the Ising like HH coupling are negligible. 
This indicates the significance of the CB admixture re- 
garding the changes of I/T2 due to strain. 
The outline for the paper is as follows. In Sec.[ll]we intro- 
duce the 8x8 k-p Hamiltonian describing states confined 
to a strained quantum dot and calculate the hybridized 
eigenstates of the HH subsystem. We find an effective 
Hamiltonian which describes the hyperfine interaction of 
the hybridized HH spin states with the surrounding nu- 
clei in Sec. |III[ In Sec. |IV| we derive the dynamics of the 
transverse HH pseudospin states and examine the effect 
of strain on the decoherence time T2 and on the hyperfine 
coupling constants. A summary can be found in Sec. W\ 
Technical details are deferred to the Appendix. 



II. HEAVY HOLE STATES IN STRAINED 
QUANTUM DOTS 

We use the 8x8 k-p Kane Hamiltonian Hk which de- 
scribes the states of bulk zincblende semiconductors in 
the lowest CB and in the HH, the LH and the SO band 
of the YJ^SHM^ We assume a flat, cylindric dot geom- 
etry which is taken into account by choosing harmonic 
confinement V^onf with lateral and vertical confinement 
lengths L and a, respectively, satisfying L ^ a. For 
detailed expressions of Hk and Vconi see the Appendix. 



Strain is added perturbatively to the sy stem by employ- 
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where the relevant block matrix elements are 
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The entries read E — ^/2P[kxexx — ikyeyy], F — Dd Tre- 



^yy ~ '^^zz), G — o/2C4[kz(£x 



'VVJ 



.) + ikyisxx - Ezz)], and J 



1/3D„(£, 
/ = Vi/2Ci[kx{e 
l/\'3 Du{£xx — Syy)- Here, P is the matrix element of 
the inter-band momentum as defined in Ref. 1251 and e^i, 
i — x,y,z, are the diagonal components of the strain 
tensor. D^ and Z)„ denote deformation potentials and 
the constant C4 is defined in Ref. UHl For simplicity, 
we restrict ourselves to a diagonal strain tensor e, since, 
due to symmetry, the shear strain components are only 
of appreciable size at the dot interfaces and negligible 
everywhere else. This assumption is valid because, due 
to their small effective mass, holes are strongly confined 
to the center of quantum dotg^S. In the vicinity of the 
F-point, the basis states in the single bands of the un- 
perturbed Hamiltonian are given by 



I* 



3,nl 



0"(r)|u±(r),±, 



(3) 



where j — CB, HH, LH, SO is the band index and ± dis- 
tinguishes between the two states of each band which are 
degenerate in the bulk system. The basis functions of i^K 
consist of s- and p-symmetric Bloch states |u~^(r)) in the 
CB and VB, respectively, and spin states |±j). The en- 
velopes are given by the three-dimensional eigenfunctions 
of the harmonic confinement potential 14onf, '/'?(r), with 



yTlx , Tly 



being a vector of the according quan- 



tum numbers. Motivated by the large energy splittings 
in quantum dots we choose nx,ny,nz G {0, 1}. We ap- 
proximately block-diagonalize the complete Hamiltonian 
H = Hk + V^onf + Hs in the HH subspace by a Schrieffer- 
Wolff transformation H = e~^He^. The eigenstates of 
the diagonal HH subsystem are determined by 



|Vl/hyb) = |^)^(l-AW)|*), 



(4) 



where A^^> is the anti-hermitian, block off-diagonal ma- 
trix describing the Schrieffer- Wolff transformation to first 
order. Explicitly, the hybridized eigenstates of the effec- 
tive 2x2 HH Hamiltonian read 



l*hyb 



(^)) = -^ E <n (^)l*J.n), 



(5) 



Table I: Values of InAs parameters we use as input for the 
8x8 Hamiltonian Hk + Konf + He . 



±, with Xlni^) 



(*I,nl^l*?JH,0>/(^.,n 



-Ehh.o) being overlap matrix elements, where H and -Ej,n 
introduce the strain dependence. i?j,n is the eigenen- 
ergy of the state Ivl/In) and M ensures proper normal 
ization. In the zero strain case we find for |\1/ 

the leading coefficients IA^b (q 1 o)(0)l = I^cb,(i,o o)(0)I 
0-11, |Alh(oi,i)(0)I = I^LH(ioi)(0)l ^ 0.097, and 



hyb(O)) 



l^so,(o,i,i)WI = l^so,(i,o,i)(0)l - 0-031. For aU con- 
figurations, Ay^Q(£) = 1. The system parameters used 
in the calculations are listed in Table H] 



III. EFFECTIVE HYPERFINE HAMILTONIAN 
OF THE HEAVY HOLE SPIN 

The hybridized HH states couple to the fcth nucleus 
by the Fermi contact interaction hi , being non- negligible 
due to the s-symmetric CB admixtures, the anisotropic 
hyperfine interaction /12 , and the coupling of the orbital 
angular momentum (0AM) to the nuclear spins ft.3 (see 
Refs. I17)29p . We derive an effective, strain dependent 
hyperfine Hamiltonian in the HH subspace by taking 
matrix elements over a single Wigner-Seitz (WS) cell: 

For the numerical evaluation of the matrix elements we 
model the WS cell as a sphere of radius one half of the 
In-In atom distance, centered in the middle of the InAs 
bond. The basis functions of i/x, |u^(r),±j), are writ- 
ten as products of 0AM eigenstates and spin states^S. 
We approximate the eigenstates of 0AM, S, P^, and 
P^ , as linear combinations of atomic eigenfunctions^^l 



+ 



],n,T 



u- (r) = a^i"^'" (r + d/2) ± Vl - a^^f^ (r - d/2) , where 
a is the electron distribution between the two atoms and 
-0"'™(r) = Rni{r)Y"'-{d, (p) are hydrogenic eigenfunctions 
with quantum numbers n, I, and m. The radial part of 
the wavefunction depends on the effective central charge 
Zcs of the nuclei where we use values for free atoms'^ ^ ^^l, 
r ± d/2 denotes the position of the hole with respect 
to the nuclei located at ±d/2 in the WS cell, where 
d = ainAs(l, 1, l)/4 is the InAs bonding vector defined by 
the lattice constant ajnAs- The bonding and anti-bonding 
character of the VB and CB are expressed by the + and — 
signs, respectively, and J^g d'^r|u^(r),±jp = 2 enforces 
normalization^. The error of this method is small and 
has been estimated in Ref. 1171 For strain distributions 



in the vicinity of the reahstic strain configuration of a 
cylindric InAs quantum dot, i.e. £xx = £yy = —0.06 and 
e^z — 0.06 (see Ref. HJ), we find an effective hyperfinc 
Hamiltonian of the form 



i^cff = {K + h')S'' + ]^{h+S- 



+ h-S+). (6) 



Here, the term proportional to hz — gh^J■BB accounts for 
the Zeeman splitting due to a magnetic field B along 
the growth direction, with g^ ~ 2 being the HH g fac- 
tor and iiB the Bohr magneton. The components of 
the Overhauser field read h^'^ — ^j. A^' {^)Ik ' where 
^^■^(e) = ^. ViAl' (e) denote the corresponding strain 
dependent hyperfine coupling constants weighted by the 
nuclear abundance Vi of each atomic species i. S is the 
pseudospin f /2 operator of the hybridized HH states and 
Ik is the nuclear spin operator of the fcth nucleus. We 
find for the effective hyperfine coupling 



Al{e) ~ v^Azie) 



funi^u)?: 



A^ie) 



E 

ij',n,n' 



z;oA±,,-,,(£)0"(rfc)*0"(rfc) 



(7) 
(8) 



where Wg is the volume occupied by a single nu- 
cleus. Az{£) and A±^j,ji{e) are the hyperfine cou- 
pling strengths and are given by Aj,ji{e) ■ I'' — 

where k, k' — ± and I*^ is the nuclear spin operator. 
In Eq. (|8|), we neglect contributions where A±jji{e) 
is more than one order of magnitude smaller than the 
leading term. We find |Az(e)| :s> ina.x\A±jj>{e)\, thus 
-ffoff is predominantly of Ising form with additional small 
pair-flip processes between nuclear and hole spin. 



IV. EFFECT OF STRAIN ON THE HEAVY 
HOLE SPIN DYNAMICS 

For a Hamiltonian of the form of Hcs, the time evolu- 
tion of the S^ (t) component and hence the decoherence 
of the HH pseudospin state is described by the Nakajima- 
Zwanzig master equatiorPSI. We obtain an algebraic form 
in the rotating frame with frequency w„ by performing 
a Laplace transform, /(s) — J^ /(i)e~'**di, Re[s] > 0, 
yielding 



S^{s + iuJn) 



(5^ 



+ I](s + iuJn) 



(9) 



Here (S'^)o = TtS~^p with density operator p and S(s) 
is the Laplace transformed memory kernel which de- 
scribes the dynamics of S^ and is derived in Refs. I10I34I 
The Zeeman splitting ujn is determined by the eigenvalue 
equation w„|n) = {bz + ft.^)|n) = {ghPsB + pAz{e)I)\n), 
where p {\p\ < 1) is the polarization of the nuclear spins 
in positive z direction, and \n) is a narrowed nuclear spin 
state^. The exact Eq. M can only be solved perturba- 
tively by expanding I](sj in powers of the flip-flop pro- 
cesses V = {h'^S^ + h^S~^)/2. This is possible since 



the energy scale of V is much smaller than the one as- 
sociated with the Ising term ^ h^ in Eq. ^. Follow- 
ing Ref. inn we expand T,{s) up to fourth order in V, 
S(s) = E(2)(s) + S(4)(s) + 0{V^), where the Zeeman 
splitting between the HH and nuclear spins forbids pro- 
cesses of odd order. E'^^Hs) and 'S^'^^s) are given explic- 



itly in Eqs. (|bI]) and (|B2| of the Appendix. E^'^^s) is 
purely real and hence leads to no decay in Eq. Q but 
to a frequency shift Aw = —Re [S(^^(s -I- iw„)] . This re- 
flects the fact that energy conservation forbids the real 
flip of the electron spin, and only virtual flips are pos- 
sible. The imaginary part of S'-'*-'(s) yields a decay, re- 
sulting in the decoherence rate I/T2 given by the relation 
I/T2 = -Im [5]('')(ia;„ + iAw - 0+)] , where 0+ is a pos- 
itive infinitesimal. T,^'^\s) and E'^^(s) are evaluated in 
the continuum limit (see Appendix). We simplify the 
calculations by averaging over the vertical dependence of 
the coupling constants A^' (e), which is possible since 
a ^ L. The frequency shift Aw can be calculated di- 
rectly, whereas the lengthy calculation of the decoher- 
ence rate I/T2 can be found in the Appendix. After a 
calculation analogously to Ref. HOI we fincpll 



I 
T2 



c+c- \A±\ 
' 4C.2 MJ 



dxxpna;] {x ^ rf)\b\{x — vi)y 



(10) 



where c± — I{I + f ) — {{ni(m ± 1))) with nuclear spin / 
and m = —I, ...,/. The brackets ((. . .)) denote averag- 
ing over all eigenvalues m of /|. 77(e) — Auj/\Az\ oc l/w„ 
and I/T2 can be evaluated numerically for any 77 < I. 
It is evident that the Ising-like form of the Hamilto- 
nian (|6]), which corresponds to \Az\ » |^±|, prolongs 
T2. The effect of non-zero strain configurations on the 
hyperfine decoherence rate is clearly visible when com- 
paring with the zero strain case. In Fig. IT] we display the 
decoherence rate I/T2 as a function of the Zeeman energy 
ujn — ghfJ-sB + pAz{e)I for an unstrained dot and dif- 
ferent internal strain configurations. The general shape 
of the decoherence rate remains unchanged for the differ- 
ent strain distributions, but the rate is shifted along the 
a;„-axis, the width of the peak is altered, and a lowering 
of the rate's maximum for asymmetric lateral strain is 
induced. The lower bound of the rates on the w^-axis is 
determined by 77 = I. The decoherence rate is shifted to 
the left with respect to the zero-strain curve for Tr e < ^ 
and shifted to the right if Tre > ^, respectively, where 
^ is a small negative number of 0(10""^). The latter re- 
lation corresponds to a dominant vertical strain tensor 
component. When the rate is shifted to the left, the 
peak becomes more pronounced, hence the sensitivity of 
I/T2 to changes in the external parameters is increased. 
The associated hybridized wavefunctions show a grad- 
ual lowering of the admixtures of all leading components 
with respect to the zero-strain case. A broadening of the 
peak occurs when the curve is shifted to the right. Here, 
the CB admixture of the hybridized wavefunction is in- 
creased while the other admixtures are lowered again. 
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Figure 1: The decoherence rate 1/72 for difFerent internal 
strain configurations as a function of tfie Zeeman energy a;„ . 

(a) For most configurations tfie maximum is sfiifted to the 
left with respect to the unstrained case except if Tre > ^. 

(b) Large ujn, i.e. magnetic fields B or polarizations p, cause 
a power-law decay of I/T2. The legend in the lower panel is 
valid for both plots and denotes the strain configurations. 



For the chosen dot geometry, L — lOnm and a = 2nni, 
the minimal coherence time at the peak of the curves is 
T2 ~ 7/iS. For large magnetic fields B or polarizations 
p, the curves decay following a power-law as evident in 
Fig. ITb). As a general result, we state that, regardless 
the strain configuration, the decoherence rate I/T2 can 
be decreased over orders of magnitude by relatively small 
changes of the external parameters. 
The strain fields in a quantum dot can be modified by 
applying additional strain, e.g. by the technique demon- 
strated in Ref. [3S1 Here, a GaAs sample containing 
InAs quantum dots is tightly glued on top of a piezo- 
electric stack, its stretching direction aligned with the 
(110) crystal axis, and a voltage is applied. So far, addi- 
tional strain of about e|| ~ 0.003 (see Ref. 155]) has been 
reached whereas the breaking point of GaAs corresponds 
to a strain of ey « 0.012 (see Ref. l37|) . We examine the 
peak of the decoherence rate located at w„.inax which is 
determined by the implicit, strain dependent equation 



dxxpnx]^ ln(a; — 77) 



-.^ 



X j ln(a; 



rj) + 7] 
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(11) 
Additional strain alters ojn^max significantly, as displayed 
in Fig. [2] and hence infiicts measurable changes on the 
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Figure 2: Shift of the peak of the decoherence rate located 
at (jjn^max whcu external strain ey up to the breaking limit 
is applied. In both panels the change is roughly linear. For 
the asymmetric lateral strain configuration (b) the relative 
change of a;„,niax is considerable larger than for the symmet- 
ric configuration (a). The legends display the internal strain 
configurations of the dots, e^.^ and e^y, at e\\ — 0. For all 
configurations we keep £zz = 0.06 constant. 



decoherence rate. The relative shift of the peak for the 
asymmetric lateral strain configuration is about 37%, 
thus larger than for the symmetric configuration where 
the shift is 22%. The strain-induced change of the deco- 
herence rate is directly connected to variances in the mag- 
nitude of the coupling strengths ^2(6) and A±jj/(£). In 
Fig. [3] we display the dependence of the absolute values 
of the dominant coupling strengths on applied strain e|| . 
We find that the relative change of A±,cb, ranging up to 
21.5%, is the largest of all whereas the relative change of 
^±,LH is only about 2%. The Ising-like coupling A^ (not 
on display) changes less than 1% and thus is negligible. 
From this we deduce that the strain-induced changes of 
the decoherence rate can be attributed mainly to the dif- 
ference in the CB admixture of the hybridized HH states. 
Hence the usually neglected contribution of the CB to the 
hole spin dephasing is of significance. 



V. SUMMARY 

In conclusion, we investigated the dynamics of hy- 
bridized HH spin states confined to self-assembled and 
hence strained semiconductor quantum dots. By taking 
into account hyperfine interaction between these states 
and the surrounding nuclei an efi^ective, strain dependent 
Hamiltonian was found, which is, for realistic strain con- 
figurations, predominantly of Ising form. The time evo- 
lution of its S^ component was derived for a narrowed 
nuclear spin state, and we have shown that the internal 
strain fields of self-assembled quantum dots affect the 
decoherence rate XjT^ significantly. For all strain con- 
figurations, I/T2 was found to be tunable over orders of 
magnitude by adjusting external parameters. Different 
strain fields were shown to cause a shift of the depen- 
dence of I/T2 on external parameters w„ and to change 
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Figure 3: The largest contributions to the hyperfine coupling 
as a function of applied external strain e\\ . The relative change 
of the CB coupling (a) is much larger than the change of the 
LH coupling (b). In the legend the internal strain configura- 
tions of the dots, e°2, and Syy, at ey — are displayed. Again 
we keep Szz = 0.06 constant for all configurations. 



272)fc2], B = \/V2P{k^ - iky), C = ;iV(2mo)73fcz(fc, - 
iky), and D = /i^/(2mo)[72(fc^ — k"^) — 2i-f3kxky], with H 
being Planck's constant, ttiq being the bare electron mass 



and 7i 



1, 2, 3, denoting the Luttinger parameters. In 



Hk, terms proportion al to C k , Bjy and B^^ were omitted 
due to their smallnesa^^I^ni^ xhe harmonic confinement 
potential 



Konf,j(r) 



"^J.-L^j.-L ,2 



-ix' + y') 1 



(A3) 
is defined by the confinement lengths L and a via 
ujj^± = h/{mj,±a^) and w^ji = h/{mj\\L'^). The effec- 
tive masses in the single bands are given by m.cB,_L/|| = 
m', mHH/LH,± = ™o/(7i T 272),tohh/lh.|| = "^o/(7i ± 
72) and msQj_/H = TOo/71. l2x2 denotes the 2x2 unit 
matrix. Table ffl in the main text lists the material pa- 
rameters of InAs we use in the calculations as input for 
the 8x8 Hamiltonian Hk + T^onf + H^ . 



the non-monotonicity of the peak. Additional applica- 
tion of external strain inflicted measurable changes upon 
l/r2 which could mostly be attributed to large alter- 
ations in the effective hyperfine coupling of the CB ad- 
mixture. This finding indicated the importance of the 
CB admixture of the hybridized HH states regarding the 
sensitivity of I/T2 on strain. 
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Appendix A: Hamiltonian - Explicit form 

In this work, we use the 8x8 Kane Hamiltonian as 
given in Ref. 1251 Appendix C, and a harmonic confine- 
ment potential. For completeness, we display here the 
general structure and relevant parts. The Kane Hamil- 
tonian is given by 
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of which the blocks 
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are relevant for the calculations, 
denoted by A = -;i^/(2mo)[(7i 



The single entries are 

h72)(fc? + fc^) + (7i- 



Appendix B: Continuum limit of the memory kernel 

For the second and forth order in a homonuclear sys- 
tem rotating with w„ we fincP^, 
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where we dropped the strain dependence of Aj^ '^ 



A^,'^(e) for readability. In both equations, c± — I {I + 
1) — {{m{m ± 1))), where / is the nuclear spin, m = 
—I,..., I and the brackets {{■■■)) indicate averaging 
over all eigenvalues m of /|. By taking the contin- 
uum limit Wo X]fe ~ /'i'^'' ^^ replace the sums by in- 
tegrals. The strain dependent frequency shift Acj = 
-Re [l](2)(s-|-iw„)] -- 10"^^eV^/w„ can be evaluated 
for Eqs. (l7| and (IS]) in a straightforward fashion. A recal- 
culation of the result of Ref. [20] gives exact shape of I/T2 
after some transformations of its original form. Starting 
from 



SW(s + iw„)~-i 
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± |2| 4± 12 



At 



' + i(Ak-AtJ' 



(B3) 



we first simplify by inserting the z averaged coupling, jus- 
tified by i ^ a and then performing the two-dimensional 
continuum limit. In the resulting two-dimensional equa- 



(A2) tion we shift to polar coordinates Xj 



i cos (fij and 



: sin ipj . After the angular integration we rescale 

i-jL. In the re- 



the radial variables by replacing rj 
suiting integral we substitute e"''" — x, rf^ = (Inx)^, 
(Iny)^ and take only the leading 



y and r|2 



term: 



from which follows 



s(4)(. 



,c+c_|A±|4 "1 "1 



A^IA, 



AxAy 



(lna;)^(ln j/)^a;j/ 
s/A,_ + i{x - y) ■ 



(B4) 



where A^ = Az{e) and A± denotes the sum of all 
A±^j.ji{e) contributing to the leading term. To calculate 
the dccoherence rate, we have to take into account that 
I/T2 — — Im [S(2cj„ + iAw — 0+)] and hence consider 



i;(^'(ia;„ + iAcj-0+) 



AulA, 



AxAy 



Jo 

2/-!^ „,\2„ 



(lna::)^(lny)^a:y 



(iAw - 0+)/A^ + i{x - y) ' 



(B5) 



Since we are only interested in the imaginary part of the 
equation, we use the following relation 



x^o C,±ix C 



(B6) 



where %, C are real numbers and V indicates that in any 
following integration of the above expression the principle 
value of the integral has to be taken. We find 



Im[ SW (ia;„+ iAw - 0+)] 



A^IA, 



da;dy(lna:;) x 



1 



-Im 



sW(z 



_c+c- \A±\ 



hzAw 

1 rl 



0^ 

AxAyxy 



(B8) 



"^^n \Az\ JO JO 

{\nxf{\nyfS{x-y-Acj/\A,\) 



(B9) 



Ax [ln.T] [ln(.T — 7])] 



c+c- \A±\ 
x{x — ri)Q{x — 77)0(1 — a; — 77), 



(BIO) 



where we use A^ — —\Az\ and replace 77 = Acj/jA^j in 
the final step. We arrive at the following expression for 
the decoherence time 



1 



c+c^ 



\A, 



^^l \A. 



Axx^nx] (a; — 77)[ln(a; — 7y)] , 

(Bll) 



(In y) y{—ii)5{x ^ y + Aco/Az), (B7) which can be evaluated numerically for any 7; < 1. 
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